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Abstract Coexistence mechanisms for species with similar ecological traits and overlapping geographic distributions 
are basic questions in ecology and evolutionary biology. Specific habitat requirements often limit distribution range 
as well as facilitate partitioning resource utilization in ecological similar species. Understanding niche segregation 
and differences in microhabitat utilization can contribute to identifying coexistence mechanisms between parapatric 
species. Feirana quadranus and F. taihangnica are two closely related frog species with parapatric geographic ranges 
and an elongated contact zone within the Qinling Mountains, which is an important watershed for East Asia. Here, 
we analysed the difference in microhabitat utilization between the two frog species and explored the key ecological 
factors that induced their microhabitat differentiation based on quadrats sampled in the contact zone. Our comparison 
of twenty environmental variables showed that both species used microhabitats with alkalescent warm water and gentle 
slope conditions. The principal component analysis indicated that climate-related variables, vegetation conditions, and 
river width were the important factors for microhabitat utilization of these species. These findings contribute to our 
understanding on the coexistence mechanisms of these two related and parapatric Asian mountain frog species. This 
study can also be helpful for identifying target habitats to conduct conservation actions and management strategies 
effectively in the face of environmental changes. 
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1. Introduction 


How to understand habitat utilization for different species 
is considered to be a crucial question in ecology and 
evolutionary biology, and also a critical component and 
the first step in species conservation and population 
management (Davidson-Watts et al., 2006; Wiens, 
1989). Various biotic and abiotic elements, such as food 
acquisition, predator avoidance, breeding, and especially 
environmental conditions (e.g. climate variability), 
contribute to habitat occupancy of species (Barber et al., 
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2008; Chailleux, 2001; Doligez and Clobert, 2002; 
Morosinotto et al., 2010). Environmental conditions are 
not homogeneous across the surface of the Earth, which 
lead to a situation that diversely combined factors are 
used by different species and demonstrate complicated 
and diversified habitat utilization by wild animals 
(Ficetola et al., 2017; Matthiopoulos et al., 2015). Given 
the influences on population dynamics, community 
structures, spatial and temporal variations of distributions, 
habitat utilization may have important ecological 
and evolutionary effects on species (Morris, 2003; 
Rosenzweig, 1981). Small changes on key environmental 
factors can alter habitat utilization of animals (Cain et al., 
2008; Luo et al., 2014). It can be reflected in both macro- 
evolutionary patterns, such as speciation and adaptive 
radiation, and micro-evolutionary strategies, such as 
characteristics of population density, persistence duration 
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and community composition (Ficetola et al., 2018; 
Morris, 2003; Rosenzweig, 1981). Divergent resource 
utilization can affect the density and persistence of 
populations, while the difference in niche occupancy may 
alter the community structure, and lead to a reduction in 
competitive interactions and promote species coexistence 
(Ficetola et al., 2018; Rosenzweig, 1981). 

Divergence in habitat utilization of congeneric species 
with overlapping ranges is a classic and important topic 
for understanding the mechanisms of speciation and 
species coexistence (Siemers and Schnitzler, 2004; Violle 
et al., 2011). Although niche conservatism indicates 
that strong competitive exclusion exists among closely 
related species and competition intensity increases with 
phylogenetic relatedness, some species do demonstrate 
similar distribution characteristics (Michalko and Pekar, 
2015; Violle et al., 2011). As proposed by Odum (1983), 
congeneric species with similar diets and foraging 
behaviours must diverge on their niches (e.g. spatial, 
trophic, or temporal niches) to facilitate coexistence 
from competitive exclusion (Davidson-Watts et al., 
2006; Michalko and Pekar, 2015). Different habitat 
utilization is considered as one of the primary driving 
forces for ecological differentiation, with varied habitat 
preferences and utilization strategies between parapatric 
species allowing them to coexist (Pita et al., 2011). 
Despite habitat utilization and preferences having been 
investigated across different scales for a various species, 
many aspects of the complicated coexistence mechanisms 
remain unclear (Cain et al., 2008; Michalko and Pekar, 
2015). 

Compared with other terrestrial vertebrates, amphibians 
have limited dispersal capability and show rather poor 
adaptive responses to environmental changes (Li et al., 
2013; Wells, 2007), though many congeneric species 
possessing parapatric ranges (Fei et al., 2009). Therefore, 
amphibians are excellent models for exploring ecological 
separation and habitat utilization differentiation (Cloyed 
and Eason, 2017; Ficetola et al., 2018). Moreover, habitat 
change is one of the most important causes of amphibian 
population declines and local extinctions (Cloyed and 
Eason, 2017; Jiang et al., 2016; Silvano and Segalla, 
2010). Understanding the differences and similarities of 
microhabitat utilization can help to grasp the ecological 
differentiation of parapatric amphibian populations, 
enrich our knowledge on their population dynamics and 
community structures, and provide important insights into 
conservation plans (Bishop et al., 2012). 

Feirana frogs in the family Dicroglossidae are 
distributed across the Qinling-Daba Mountains (QDM) 
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(Fei et al., 2009; Hu and Jiang, 2018; Wang et al., 2009). 
Based on phylogeographic analyses, the ancestors of Æ 
quadranus (swelled-vented frog) radiated in the QDM and 
then expanded to the Longmen-Micang-Daba Mountains, 
while F taihangnica (Taihangshan swelled-vented frog) is 
native to the Qinling Mountains (Wang ег а!., 2009, 2012, 
2013). A substantially elongated contact zone between the 
two species is located in the central Qinling Mountains 
(Wang et al., 2009; Yang et al., 2011). It is suggested that 
weak environmental tolerance restricted the dispersal of 
Е quadranus further north, and interspecific competition 
was likely preventing the southward expansion of 
F. taihangnica (Hu and Jiang, 2018). Additionally, 
these parapatric species have similar natural histories 
and closely resemble each other ecologically and 
phenotypically, except for one significantly different 
trait of breeding males (i.e. Е taihangnica possesses a 
multitude of tiny granules on the swollen skin of the anal 
area but Е quadranus does not) (Fei et al., 2009). Thus, 
the two frog species present a remarkable opportunity to 
study the microhabitat segregation of congeneric species. 
In this study, we investigated the differences and 
similarities in microhabitat utilization between Е 
quadranus and F. taihangnica. Specifically, our aims 
were to examine: (1) what different characteristics exist 
on microhabitat utilization between these two species, 
(2) how they share homogeneous resources in local 
assemblages, and (3) what ecological factors influence 
their microhabitat utilization. 


2. Materials and Methods 


2.1. Study Area The study area is located in the Qinling 
Mountains (33°16’-34°0S'N, 107°20'-109°04’E; Figure 1), 
which forms the transition area of subtropical and warm 
temperate zones in China and is an important watershed 
for East Asia (Zhang et al., 2017). The north slope of the 
Qingling Mountains has a warm temperate and semi- 
humid climate because of the cooling dry air from the 
northwest, while the south slope has a subtropical and 
humid climate under the effects of the humid air from 
the southeast (Zhang et al., 2017). There are mainly three 
vegetation types along the altitudinal gradients (from low 
to high): subtropical evergreen and deciduous broadleaf 
forests, temperate deciduous broadleaf and subalpine 
needle leaf forests, and subalpine scrub meadow (Loucks 
et al., 2003). Our study was carried out within the region 
between the Wei River and the Hanshui River in Shaanxi 
and Gansu provinces, which locates in the contact zone 
of Е quadranus and Е taihangnica (Figure 1) (Hu апа 
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Jiang, 2018; Wang et al., 2009, 2013). 


2.2. Data Collection Based on drainage conditions 
of the study area, the known distributions of Feirana 
species, and the geographical locations of main rivers 
in the QDM, seven main tributaries of the Wei River 
and the Hanshui River were chosen to conduct transect 
population investigations: the Fengyu, Shitou, Qianyou, 
Jiaoxi, Xushui, Hongyan and Yanchuan River (Figure 1). 
They contain different characteristics to be representative 
of most rivers in the overlapping region. Our field 
investigations were conducted from August to September 
in 2010. We searched for the individuals of target frog 
species along the transects, and 1 т х 1 т quadrats were 
set as sampling points to measure environment variables, 
with the individual location as the quadrat center. Then, 


100°E 


105°E 


— Toa Meer | 


Sf 


110°E 


referring to the previous study (Buckley and Jetz, 2007) 
and our field experience on target species, 20 environment 
variables, including air temperature, water temperature, 
humidity, water depth, pH, slope of riverbank, riparian 
canopy, distance to water, distance to road, current 
velocity, river gradient, river width, river width in 
flood season, sandy component, vegetation coverage, 
position on slope, and vegetation type of each quadrat 
were measured and recorded (Table 1). These variables 
may impose constraints on physiology, distribution, and 
survival of amphibians (Wells, 2007). We measured 
48 quadrats for Е quadranus and 62 quadrats for F. 
taihangnica, respectively. Since there was only one 
quadrat that possessed both of the two species, we did not 
consider niche separation index in this research. 
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Figure 1 Geographical location of the study area and sampling quadrats for Feirana quadranus and F. taihangnica in the Qinling Mountains. 
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Table 1 Environment factors compiled to depict the microhabitat utilization of Feirana quadranus and Е taihangnica in the Qinling 


Mountains. 


Environment factors 


Descriptions 


Air temperature (°C) 
Water temperature (°C) 
Humidity (%) 

Water depth (cm) 

pH 


Slope of riverbank (°) 


Distance to water (m) 
Distance to road (m) 


Riparian canopy (%) 


Current velocity (m/s) 


Position on slope 
River gradient (°) 


River width (m) 


River width in flood season (m) 


Sandy component (%) 


Vegetation type 


Air temperature of the centre of each 1 m х 1 m quadrat 

Water temperature of the centre of each 1 m х 1 m quadrat 

Humidity of the centre of each 1 т х 1 т quadrat 

The average value of the deepest and shallowest depth of water in the quadrat 


The average value of three measurements of the centre of each 1 m х 1 m quadrat by using a pH meter 


The slope of riverbank in each | m x 1 m quadrat by using the DQS-2A theodolite, L for the left bank and R for the 
right bank 


The distance of each quadrat centre from the nearest deep water (water depth > 50cm) 
The distance of each quadrat centre from the nearest road (artificial road, road width > 1m) 


Stand in the centre of the quadrat, estimate the riparian canopy above the water 


The time that a self-made foam drifting for 1 m in the water, and the average value of three measurements was used 
to represent the water velocity 


The position of each quadrat on the slope, and it was recorded with three categories: upper (1), middle (2) and lower 
(3) 
The slope of the river in the quadrat centre by using the DQS-2A theodolite 


The width of the nearest stream stretch from the centre of quadrat 
The width of the nearest stream stretch in the flood season from the centre of quadrat 
The proportion of sand in the sediment of a river in the quadrat 


The types of vegetation on the left and right river bank, respectively. The vegetation type in each quadrat was 
determined by primary plant species, and recorded with six categories: arbour (1), shrub (2), herbage (3), arbour and 


herbage (4), arbour and shrub (5), shrub and herbage (6); L for the left bank and R for the right bank 


Vegetation coverage (%) 


2.3. Statistical Analyses We tested normality of each 
continuous variable using a Kolmogorov-Smirnov test, 
and transformed the percentage data by arcsine/square- 
roots and made In-transformations for variables that 
were non-normally distributed. Air temperature, slope 
of riverbank, distance to water, and river width were In- 
transformed. If the variable data met the assumptions of 
normal distribution and equal variance, we performed an 
independent-samples 7 test to compare its means between 
Е quadranus and Е taihangnica; otherwise, we used a 
Mann-Whitney U test. Additionally, we conducted Chi- 
square tests to assess whether there were significant 
differences of the categorical variables (i.e. position on 
slope and vegetation type) between species. 

To reduce the potential collinearities among variables, 
we conducted Spearman’s correlation tests and excluded 
humidity data in the later analyses (highly correlated 
with air temperature, |r| > 0.70; Supplementary 
Materials, Table S1). Then, we conducted a principal 
component analysis (PCA) to analyse the effect of each 


The percentage of aboveground vegetation cover of the quadrat; L for the left bank and R for the right bank 


ecological factor on microhabitat utilization of the frogs. 
Principle components (PCs) that had eigenvalues > 1.0 
were extracted for subsequent analyses. We used the 
independent-samples 7 tests (for normal data) or the 
Mann-Whitney U tests (for abnormal data) to detect the 
interspecific differences of the PCs. All analyses were 
conducted using SPSS 16.0 (SPSS Inc., Chicago, IL). 


3. Results 


In general, both of the frog species used flat habitats (river 
gradients: Е quadranus/F. taihangnica, respectively, 
7.09 + 4.05°/6.18 + 5.18°; mean + SD; the same below), 
with relatively warm (water temperature: 18.72 + 4.32 
°C/17.36 + 3.72 °C) and alkalescent water (pH: 8.36 
+ 0.63/8.13 + 0.45; Table 2). Е quadranus occupied 
microhabitats with significant higher humidity, steeper 


slope of the left riverbank, and faster current velocity 
than those of Е taihangnica (T-tests, all Р < 0.05; Table 
2). More sandy component and less vegetation coverage 
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Table 2 Results of the Mann-Whitney И tests and Chi-square tests, showing the differences of the ecological factors affecting microhabitat 


utilization of Feirana quadranus and Е taihangnica. 


Е quadranus 


Е taihangnica 


Ecological factors t/Z P 
Mean SD Mean SD 
Air temperature 3.251 0.198 3.262 0.202 –0.310 0.757 
Water temperature 18.715 4.321 17.356 3.723 –1.396' 0.163 
Humidity 1.16 0.196 1.042 0.193 3.147 0.002" 
Water depth 2.791 0.538 2.894 0.542 —1.067 0.286 
PH 8.352 0.631 8.134 0.447 –1.396' 0.163 
Slope of river bank (L) 3.575 0.664 3.055 0.838 3.529 0.001" 
Slope of river bank (R) 3.608 0.602 3.441 0.778 1.232 0.221 
Riparian canopy 0.391 0.394 0.492 0.456 –1.101' 0.271 
Distance to water 1.784 1.217 1.406 1.16 1.655 0.101 
Distance to road 30.604 32.754 24.632 23.043 ~0.637' 0.524 
Current velocity 0.414 0.23 0.25 0.245 3.594 < 0.001" 
River gradient 7.09 4.052 6.181 5.177 –1.862' 0.063 
River width 1.675 0.802 1.515 0.78 1.049 0.297 
River width in flood season 21.4 20.293 21573 18.286 –0.263' 0.792 
Sandy component 0.59 0.236 0.451 0.375 -2.293 0.022" 
Vegetation coverage (L) 0.902 0.315 0.986 0.304 –1.847' 0.065 
Vegetation coverage (R) 0.961 0.326 1.067 0.261 –2.153! 0.031" 
Position on slope 2.607" 0.272 
Vegetation type (L) 6.954" 0.224 
Vegetation type (R) 7.757 0.17 


Note: L, the left bank; В, the right bank; ° Independent-samples T tests; | Mann-Whitney U tests; * Chi-square tests; * Significant difference 


between species. 


of the right riverbank were found for F quadranus than 
Е taihangnica (Mann-Whitney U tests, both P < 0.05; 
Table 2). No obvious interspecific difference existed 
in the position on slope or the vegetation type on both 
riverbanks (Table 2). 

Our PCA extracted seven PCs which explained 66% 
of the total environmental variance, with the Ist and 2nd 
PCs explaining 19.6% and 11.2%, respectively (Figure 
2; Table 3). According to the component matrix, water 
temperature, air temperature, pH, humidity, riparian 
canopy, river width, river width in flood season, and 
vegetation coverage were the key factors that influenced 
microhabitat utilization (Table 52). The first five PCs 
contained the information in majority of the variables 
(55%) and were compared for the interspecific differences 
(Table S1). Limited differences were displayed in all the 
first five PCs between F quadranus and Е taihangnica 
(Table 3). 


4. Discussion 


Parapatric species usually occupy divergent niches, 


44 OF quadranus А 
А F. taihangnica 

24 
= А 
x 
а о 
= 
= 
5' о 
Е о 

о 
224 о 
2 о 
5 о 
А 
44 А 


РСІ (19.6%) 


Figure 2 Microhabitat utilization of Feirana quadranus and 
F. taihangnica based on the first two PC axes. We displayed 10 
variables with component loadings > |0.50|, including WT (water 
temperature), AT (air temperature), pH, VCR (vegetation coverage 
of the right bank), RW (river width), VCL (vegetation coverage 
of the left bank), RC (riparian canopy), POS (position on slope), 
SRBR (slope of the right river bank) and SRBL (slope of the left 
river bank). 
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Table 3 Results of the principle component analysis (PCA) for environmental factors and the independent-samples 7 tests of PCs between 
Feirana quadranus and Е taihangnica, reporting the main contributors and no significant inter-species difference of the PCs. We only present 
variables with component loadings > |0.50| and the first five PCs which contained the information in most of the variables. 


Components Component loadings Explained P 
variance 
PCI water temperature (0.86), air temperature (0.73), pH (0.67), vegetation coverage (L: –0.60; R: –0.62), river 19.56% 0.064 
width (0.60), riparian canopy (—0.59) 
PC2 position on slope (0.58), slope of river bank (L: –0.56; R: –0.57) 11.19% 0.103 
PC3 vegetation type (L: 0.59) 8.90% 0.263 
PC4 river gradient (—0.56), sandy component (0.52) 7.74% 0.128 
PCS water depth (0.70) 7.26 0.087 


Note: L, the left bank; R, the right bank. 


such as segregated microhabitat, to sustain species 
coexistence over a long term (Morris, 1996; Pita et al., 
2011). Our results suggested differences in microhabitat 
utilization for the two parapatric Asian mountain 
frog species. In comparison with F. taihangnica, F. 
quadranus preferred humid microhabitats that were 
on steeper riverbanks, with less vegetation coverage, 
faster current velocities, and more sand component on 
riverbed. Since differences in habitat preference among 
individuals can lead to disparities on gene frequency 
within species, habitat complexity plays a major role 
in speciation process of congeneric species (Graham 
et al., 2004). Phylogeographic analyses showed that F 
quadranus and Е taihangnica were diverged 4.6 million 
years ago (Wang et al., 2009). Based on the deduced 
historical expansion, F. taihangnica is considered to 
be a native species while Е quadranus is probably ап 
external competitor within their overlapping ranges in the 
central Qinling Mountains (Hu and Jiang, 2018; Wang 
et al., 2012). Specific habitat requirements often limit 
distributions of animals as well as facilitate isolation of 
ecological requirements (Siemers and Schnitzler, 2004). 
It is suggested that niche differentiation can lead species 
to complement one another in mixture and, hence, result 
in better utilization of available resources (Odum, 1983). 
The difference of microhabitat utilization between the 
two frog species might demonstrate a result of variable 
individual preferences on environmental conditions after 
a long period of evolution. 

Since the two congeneric frog species originate from 
different areas and have some distinct phenotypic features 
(Fei et al., 2009; Wang, et al., 2009), it is expected that 
they can occupy different microhabitats in the contact 
zone to reduce interspecific competition and maintain 


species coexistence. However, despite slight differences 
on microhabitat utilization between the two frog species, 
interspecific divergence on ecological factors generally 
influencing microhabitat utilization was lacking in this 
study. Ecological requirements of a species often reflect 
the balance between the effects of diversified intraspecific 
competition and constrained interspecific competition, 
which can not only promote coexistence of different 
species but also hinder habitat types the animal use 
(Cadotte and Tucker, 2017; Mayfield and Levine, 2010). 
Additionally, the niche filtering theory indicates that, in 
a given environment, the profitable traits are selected by 
local conditions and some species could be filtered out 
because they cannot cope with local elements (Michalko 
and Pekar, 2015). It might then lead to a coexistence of 
species with similar ecological needs and characteristics. 
This means that competition can sometimes eliminate 
taxa with distantly genetic relationships or differently 
ecological traits (Mayfield and Levine, 2010). The limited 
segregation of microhabitat utilization in this study 
supported the expectation that closely related species 
might have similar ecological niches due to the suites of 
environments (Hu еѓ al., 2016; Warren et al., 2014). 
Leibold and Mcpeek (2006) suggested that weak 
interspecific competitions would not make parapatric 
species to narrow their current niches. This deduction is 
partly supported by the absence of extensive microhabitat 
segregation between F quadranus and Е taihangnica, 
which, might indicate that the available microhabitat 
resources have not been completely occupied by these 
frogs in the Qinling Mountains. On the other hand, 
the phenomenon of relatively large variations in some 
ecological factors (e.g. river gradients, pH, and water 
temperature) may lead to a speculation that the frogs 
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are not sensitive to those environmental variables. 
Microhabitat utilizations investigated by this study 
are indeed only a part of the realized niches of these 
frogs. Even though no obvious ecological microhabitat 
divergence was detected, niche divergence at other levels, 
such as in their prey selection or reproductive resources, 
may be at play instead. Remarkably, considering the 
amount of environmental variables used in this study, 
the number of sampling quadrats was relatively small. 
Different types of data for the environmental factors, 
such as temperature-related variables (continuous data), 
vegetation conditions (categorical data), and river 
conditions (percentage data), were included, which could 
lead to a low interpretation of the PCA and, thus, might 
cause potential uncertainties of the results. Therefore, we 
suggest that more extensive surveys should be carried out 
for understanding microhabitat utilization of these frogs. 

Our study indicated that discrepant preferences on 
environmental conditions are likely to cause microhabitat 
segregation between congeneric and parapatric species. 
Several mechanisms of species coexistence can interact, 
for example, niche partitioning on one dimension may 
facilitate niche filtering on another (Michalko and Pekar, 
2015). When availabilities are high for some types of 
resources, specific preference rather than competition 
may play roles in microhabitat utilization. However, 
due to the limited resources, filtering effectiveness on 
the environmental factors could obviously influence 
coexistence of ecologically similar species (Mayfield and 
Levine, 2010). 

Amphibians have highly water-permeable skins and are 
sensitive to variations on ambient environments (Wells, 
2007). Environmental changes can pose intense impacts 
on amphibian species. For example, sharp changes in 
temperature can affect their survival and reproductive 
processes (Benard, 2015). Studying microhabitat 
segregation between parapatric amphibian species can 
be helpful for understanding their ecological niche 
separations, identifying target habitats, and then providing 
scientific support to establish efficient management 
policies and conservation plans for those vulnerable 
species (Bishop et al., 2012). In this study, microhabitat 
utilizations for both of the species were highly depended 
on vegetation and river conditions. This suggested that 
keeping relatively high vegetation coverage and avoiding 
encroachment on natural rivers should be effective to 
protect them. Nevertheless, interspecific comparisons 
on other ecological traits (e.g. dietary, predator, and 
reproductive strategy) should be carried out in the future 
to further understand the coexistence mechanisms of 


these two parapatric species and identify the protection 
keystones. 
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Table S2 Results of the principle component analysis (PCA), explaining variations, eigenvalues, and factor loadings for the seven principal 
components with eigenvalues > 1. 


Component 
1 2 3 4 5 6 7 
Water temperature 0.855 0.208 0.108 —0.131 0.101 —0.187 –0.039 
Air temperature 0.729 0.389 0.138 –0.137 0.127 –0.220 –0.006 
РН 0.671 –0.066 0.234 –0.053 0.043 –0.191 –0.336 
Vegetation coverage (R) 0.618 0.261 0.026 0.018 0.224 0.487 0.073 
River width 0.600 0.066 –0.259 0.051 0.213 0.032 0.251 
Vegetation coverage (L) –0.596 0.229 0.198 0.092 0.240 –0.167 –0.239 
Riparian canopy –0.589 0.136 0.140 –0.174 0.140 0.328 –0.144 
River width in flood season 0.485 0.268 –0.338 0.101 0.189 0.438 –0.219 
Position on slope 0.311 0.580 0.334 –0.003 0.044 0.217 0.263 
Slope of river bank (R) 0.106 –0.573 0.328 0.107 0.304 0.247 0.149 
Slope of river bank (L) 0.258 –0.557 0.280 0.120 –0.061 0.301 0.113 
Current velocity 0.129 –0.484 0.324 –0.351 –0.211 —0.052 –0.149 
Distance to road –0.245 0.471 0.297 –0.095 0.185 0.278 0.471 
Vegetation type (L) –0.093 0.095 0.590 0.540 –0.153 0.117 –0.075 
River gradient 0.025 –0.190 0.408 –0.557 0.004 –0.252 0.371 
Sandy component 0.088 –0.162 –0.176 0.519 –0.313 –0.288 0.429 
Distance to water –0.184 0.184 0.025 –0.496 –0.344 0.187 –0.107 
Water depth 0.002 –0.181 0.297 0.184 0.695 –0.143 –0.166 
Vegetation type (R) 0.144 0.378 0.489 0.267 –0.493 –0.066 –0.202 
Eigenvalue 3.717 2.127 1.691 1.471 1.379 1.185 1.064 
Cumulative variance (%) 19.563 30.755 39.657 47.397 54.657 60.893 66.491 


Note: L, the left bank; R, the right bank. 


